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APPROXIMATE REPRESENTATION 
Br W. E. Byebly 

1. Introduction. That the Least Square criterion employed in this 
paper could be made to give the received coefficients in a Fourier or in a Zonal 
Harmonic development appears to have been first noticed by Plarr who com- 
municated his discovery through Bertrand to the Comptes Rendus in 1857.* 

Later Toepler f and Gram % who were apparently unacquainted with 
Plarr's paper used practically the same method as Plarr ; and these methods 
have been employed to some extent by recent writers on Integral Equations. § 

So far as I have been able to discover these writers have usually confined 
their attention to the one-dimensJonal case, although some of them have men- 
tioned that the methods they employed could be extended to the case of two- 
dimensional and three-dimensional problems. 

That the methods here discussed furnish a very convenient and powerful 
tool in dealing with problems in Mathematical Physics which require the use 
of the so-called " Harmonic Analysis " is not, however, generally known, and I 
have found them so useful in my own work that it seemed to me worth while 
to communicate them to my fellow workers. 

2. The Mean Square Criterion. Suppose that/(x) and F{x) are 
two given functions and that we wish to regard the second as an approximate 
representation of the first over a specified range of values for x, i.e., 
Xq < X < Xx ; is there a convenient criterion by which we can judge the excel- 
lence of the approximation ? 

When such a representation is used it is ordinarily desirable that the ab- 
solute value of the error E =f{x) — F{x), which of course is a function of x, 
should be small for all values of x considered. Of course in that case 
E^ = [/(x) — F(x)y will be small and [jE'^], the mean value of E^ over the 
range in question, will be small. 

* Comptes Jtendus, vol. 44, p. 985. 

t Wiener Anzeige, vol. 13, (1876), p. 205. 

t Crelle, vol. 94, (1883), p. 41. 

§ M. B6cher, An Introduction to the Study of Integral Equations, pp. 52-60. 

^128) 



APPROXIMATE REPRESENTATION 129 

This [S"^], the Mean Square of the Error, that is, the mean value of the 
square of the error j{x) — F{x), over the i-ange in question, we shall use as 
the measure of the excellence of our approximation. The smaller [JS^], the 
better we shall call our representation. 

Of course, since the mean of a set of positive quantities cannot be zero 
unless every one of the quantities is separately zero, our representation is 
exact if [^*] = 0, and in that case /(a;) and jF'(x)are equal for all values of a; 
between x^ and x^. 

It is easily seen that 

[j^] = ^-1^ £ [/(x) - F{x)Y dx. [I] 

3. Determination of the "Best" Coeflacients in an Approx- 
imate Representation of Specified Form. Let f{x) and ^i(x) be 
given functions* of x and let us try to determine a constant coefBcient a^ so 
that tti ^i(x) shall represent f(x) over the range from x^ to Xi as well as pos- 
sible. 

The problem is a very simple one in maxima and minima. 

[^^] ^ -L^J\f(x) - a,4>,{x)Y'dx, 
regarded as a function of a^, is to be made as small as possible. Write 
^P. ^ - -^ riAx) - a,U^)] U^) dx = 0; 

«"l Xj — Xf) jx^ 

O 
we get ai = -^, [II] 

where -^i = [<f>i{^)Ydx 

and Ci = / f(x)^i(x)dx; 

then r^^l = ^ \^ r r\f{x)fdx - 2aiCi + a\A^ . 

* Throughout this paper we shall suppose that the functions employed are of a " re- 
spectability " above suspicion; fiuite, continuous, single-valued, integrable, and difFerentiable 
over the range or throughout the region considered. It is true that in many cases some of these 
limitations may be dispensed with, "but that is another story." 
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But from [II] 2aiOi = 2a\Ai. 

Hence 

4. Having found a^ let us now determine a^ so that tlie addition of tiie 
term 0,2 4>i{x) , ^^(x) being another given function, slmll improve our repre- 
sentation as much as possible. To do this we have only to determine Oj so 
that «.j (f>i{x) shall represent /(«) — «i <l>i{x) as well as possible, and we get, 
by §3, 

f(^) (f>i{'<^)dx - aj / 4>i(x) 4>2(x)dx 

J\<f>2{x)ydx 

or, if 

Ai= (i<^i{x)fdx, ^2= f[<f>.,(x)y dx, 0,= fJix) <j>(x) dx, 

Jx, Jx„ Jx, 

fX, /-Xi 

0.i= f(x) (}>2{x)dx, and ^1,2= N>i(x)^.,{x) dx, 

Jx„ Jxo 



Oo = 



then 



«2 = 2 ' [IV] 



and [^-^j = j-A^^ [jr[/(x)]^cix - a\A, - al^]- [V] 

5. Let us go back and determine the coelBcients aj and a^ ab initio so 
that Ui 4>\{^) + *2 'A2(^") shall represent /(as) as well as possible over the 
range x^ to x^. 

Here we have a problem in maxima and minima of a function of two inde- 
pendent variables. 



write 



\eA = -i— rif(x) - a,<}>,(x) - a2<}>,{x)Ydx; 
L J Xj — X0 jx^ 



8[jE^] _ -2 
9ai 



—2 r^i 

*1 ~ ''^0 Jxn 

J-— J- = / [/(x) - ai^i(x) - a3<^2(x)] ^2(x) dx = 0. 

oa^ Xi — Xo jx„ 
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To get Oj and Oj we have to solve the linear equations 

A"! + -Bifia^ = Ci, [VI] 

where the notation is the same as in §4. 

From [VI] ia-^Gj, = 2alAi + ia^a^Bi^^, 

Hence 

\.^^~\ ^ sT^o [jy^''^^'' '^"' ~ "'^' ~ "^'^^ ~ ^"i'"^^!-^] • t VII] 

6. Example 1. (a) To determine % so that ajX shall be the best ap- 
proximate representation of sin x from a; = to x = tt. 

Here 

Ci = X sin xdx — tt, Ai — a^dx = 

Jo Jo 



3 ' 



ai = ~ = -^ = 0.304 ; 



and sin x = 0.304x approximately if < x < tt. 

As our test of the excellence of our approximation we have, since 



I 



sin* xdx = — , 



It is interesting to note that should we take for sin x the first term of the 
familiar power series 

x' x" 

'''"31+ 51 
we should have over the same range 

M=t'-I= '■'««• 
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and for 



(6) To determine a<i so that the addition of the term a^v? shall improve 
our representation as much as possible. 

Here A^^ and Gx have the same values as in (a) . 

A = y. Ca = TT* - 67r, and B^^^ = — ; 

sin X = 0.304a; - O.OlSa;^ 

M=i[5-T''!-T''5]=''-°^'"- 

(c) To determine a^ and a^ so that ajS; + a^x? shall represent sin ic as 
well as possible over the range < a; < tt. Here A^, A^, Ci, C^, and Bi^ 
have the same values as in (6) . We have 

Y^i + -5-^2 = '^' -5-^1 + y-^2 = TT* — 67r. 

Hence 

«.= ^(|i- = 0.8576, 

and 

sin a; = 0.8576a; - 0.09348x« 

approximately if < x < tt. This last approximation is excellent since 

Should we use as our approximation 

sm X = X — ■— f 

O 4 

we should have 

[B^} = 0.4008 
ver the same range. 
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7. Example 2. (a) To determine Oi so that ay sin x shall represent 
X over the range from to tt as well as possible. 
Here 

Ai= I sin^ X dx = - , Oi= I x s'm x dx — tt, 

<.. = |! = 2, 

and 

X = 2 sin X 
approximately (0 < x < tt) . 
As our test of the excellence of the approximation we have 

(6) To determine a^ so that a^ sin 3x shall improve the value just ob- 
tained as much as possible. 
Here 



so that 



/" TT 

sin^ Sx dx = — , 

Jf" TT 

f X sin 3x dx = — , 
^ 

Bi2= /sin X sin 3x dx = 0, 



"■^ = 3' 



here r„„-l ITtt* ^ 2 1 tt^ 2 ^ ^, 

(c) To determine Oj and a.^ so that Ui sin x + a2 sin 3x shall represent 
X as well as possible over the range < x < tt. We have 

TT TT TT 

-ai = -TT, 2"*= 3' 

so that 2 

ai = 2, and «2 = o ; 



M=Ut-^'-H ='■»'■ 
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8. Example 3. To determine a^ and a^ so that aj sin a: + a^ sin'ac shall 

IT 

r 

Ai= sin* X dx = — , 
Jo 4 



represent sin 3x as well as possible over the range < x < -x. 



/""/* Sir 

^2 = / sin» a; cZx = ^, 

Bii= / sm a; sin' a; aa; = r^, 

/•it/: 



rt/i 

sin 3x sin a; dx = 0, 



/■t/2 ,r 

0^= I sin 3x sin' x c?x = — ttt' 



TT 37r ^ Stt 57r -jr 



Hence «! = 3, 

a^ = —4, 

and the representation is exact, and 

sin 3x = 3 sin X — 4 sin* x 

TT 

for all values of x between and - . 

9. The method of §5 is easily generalized. Let it be required to deter- 
mine «!, a^, •••,«„ so that 

k= n 
k=:l 

shall represent /(x) over the range Xo < x < Xj as well as possible. 

We find that the coefficients can be obtained by solving the n linear 
equations 






2) A.z«* + ^zflz + 2 ^*.'«* = ^' ' » (^ = 1 . 2, 3, • . ■ n. )• [VIII] 
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( r j=" fan— 1 l=n \ 

/'a, 

where -4^ = I l4>,f(x)Y dx, 

-Si,j= /^(a;) <^,(a;) dx, 

and Ci = / y"(a;) ^^(a;) (Zx. 

If we desire to correct our representation as much as possible by the ad- 
dition of another term a„^i <^„ + i (a;), 

Cn41 — _^ «* -Si,,i+i 






••n+l 

and 

t-i 1 C r*i *=n+l *=n— 1 /=» ) 

or [^^] = [^]-A+i<+i- [XI] 

10. Orthogonal Functions. A pair of functions the integral of 
whose product taken over tlie range x^ < x < x^ is zero we shall call orthogonal 
over the range. 

That is if 

/ 4>m{x) ^ni:^') dx = 0, 

(f>,a(x) and (f>n(x) are orthogonal over the range x^ <x < Xi. A set of functions 
such that every pair is orthogonal over our range we shall call a set of ortho- 
gonal/unctions over the range. 

If in §9 the functions <f>i(x), <^2(^)> • • • ^n+i(^) ^^^ orthogonal our 
results simplify surprisingly. 

Since J^k,i= / 4>k(^) 4'i(^) dx= 
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[VIII] gives at once 

A = ^; ; [XII] 

[IX] gives 

and [X] gives 

«„ + i = -X^, [XIV] 

while [XI], 

remains unchanged. 

We see then that if we are trying to represent/ (x) as well as possible by 
the aid of a sot of functions which are orthogonal over the range considered 
the coefficient a/, of a term a^. 4>k(^') which would be the best one if we used 
that term alone is the best one if that tei*m is used in combination with any 
others of the set, and is the best one if that term is to be added as a correction 
to terms of the set whose best coefficients had been previously determined. 
Moreover since [^^] is necessarily positive no matter how many terms are 
used in the represeution it is clear from formula (XIII) that the addition of any 
term with its appropriate coefficient* if the term can form an orthogonal pair 
with every term already present diminishes the Mean Square of the Error and 
improves the representation. 

11. Corollaries. A few easy corollaries are : 

(a) Iif{x) is capable of being exactly expressed in terms of the assumed 
functions for all values of x between Xq and x^ formula [VIII] §9 applied to 
the range aso < a; < Xi or to any portion of that range will give the correct 
coefficients. If it has not been foreseen that exact expression is possible for- 
mula [IX] will establish the fact by showing that {_£!^] = 0. (Compare §8, 
Ex, 3). 

(6) If it is known that for all values of a; in a given range f(x) can be 
expressed as a series of terms of the assumed form, and all the functions 

*It must be kept in mind, however, tliat the appropriate coefficient Ct/At may happen to 
be zero, in which case of course the improvement is null. 
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(j>ii:{x) employed are orthogonal over the range, formula [XII] §10 will give 
the correct coefficients of the series. 

(c) No development oi f{x) in a series whose terms are orthogonal 
over the range over which the development is to hold good can be correct 
unless the set of forms <^a(x) employed is a complete orthogonal set; that is 
unless all forms are present which might be inserted without preventing the 
whole set employed from being orthogonal.* 

{d) Since in computing the value of a function by the aid of a series 
into which it has been developed one can use only a limited number of terms, 
the check given by the value of [^] is as valuable in practical applications in 
Physics when a valid series is employed for f(x) as when an approximate 
representation is used. 

(e) If in all that has gone before we replace function hy point-function, 
rangehy region (one-, two-, or three-dimensional) , and integral over the range 
by integral throughout the region, none of our formulas nor of our results will 
be changed in any essential particular. 

12. Trigonometric Series. Since sines and cosines of whole multi- 
ples of X form a set of functions orthogonal over the range — ir < x < ir the 
coefficients of the terms of Fourier type that best represent f{x) over that 
range are given by formula [XII] §10. 

But these are the familiar Fourier Coefficients. 

Let/(.^) = ^hy,+ bi cos x + b^ cos 2a; + • • • + % sin a; + a.^ sin 2x + • • • 

1 f" 
Then b„ = — f(x) cos mx dx, 

1 f' 
and a^ = - / f(x) sin mx dx. 

[^] = ^ j /;/(=«)]' dx--^x «^ + ^'«) I • 

Of course if /(x) is developable in a Fourier's Series our coefficients «,„ and 
b^ are the correct coefficients (v. corollary (b) §11). 

13. Let us obtain and test a few Trigonometric expressions for x over 
the range < x < tt. 

♦ Apparent exceptions to this corollary are due to the fact that the appropriate coefficient 
C'k/At for a term of proposed form may happen to be zero. 
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Here, since when m and n are integers / sin tnx sin nx dx = 0, and 

/cos mx cos 7ix dx = 0, while / sin mx cos nx dx is not generally zero, we 
Jo 

can deal conveniently with a sine expression 

X = ^S^ a^ sin mx, 

or a cosine expression 

X = ^S*' b^ cos mx. 



These are 



^rsin a; sm 2a; sin 3a; sin 4a; ~i 

TT 4r cos 3a; cos 5a; "i 

a;=^--|_cosx+-^5-+-^+ ...J. (2) 

For (1), m = l[^- 2- (i + ^+ i+ • • •)] ='' 



since 



1 1 1 


TT^ 


P "*" 2* "*" 32 "^ ' ' 


■ ~ 6 



(v. Chrystal's Algebra vol. II, page 243), and as is well known (1) is a valid 
development of x over the range < a; < tt. 
For (2) 

111 IT* 

since _ + _ + _+... = _ 

(v. Chrystal's Algebra, vol. Ii, page 243), whence 

ill - ^ 

and as is well known (2) is a valid development of x. 

' sin mx cos nx dx = 



P 
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if m + n is even, that is if m and n are both odd or both even. Sines and co- 
sines of odd multiples of x form, then, a set of functions orthogonal oyer the 
range < x < tt, as do sines and cosines of even multiples of x. We can get 

./since sin 3a; sin 5x \ 4 / cos 3x cos 5x \ ,^^ 

x = 2(-^ + _^ + -^+...)--(cosx + -^ + -^ + ...),(3) 

and x = |-2(?ilL!^ + !llL!^ + ^-i^+...). (4). 

For (3), 
^„„lr7r3„/l 1 \8/l 1 1 Mn 



since 

For (4) 

r^,^lr7r3 7r3„/l ] 1 \-ln 

and (3) and (4) are valid developments of x from x = to x = tt, which are 
not usually given in the texts, 

14. Spherical Harmonics. Our method gives the familiar coeffi- 
cients in the case of the other Harmonic developments. 

(a) Zonal Harmonics- For instance since 

/P„(x) P„(x) dx = 

if m and n are integers, the Zonal Harmonics Pct(x) form a set of functions 
orthogonal over the range — 1 < x < 1. If we wish to express f{x) as 





1 

1* 


+ 


1 
32 


+ 


1 
5'^ 


+ 


ttS 




TT^ 






r 


1 


3 


~~ 


4 


~~ 


2.. 


¥ 



^ OJm Pm {^h 



and 



y f{x) p,„{x)dx 
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The Zonal Harmonics P^C^os 6) are orthogonal over the surface of the 
unit sphere. (Byerly's Fourier's Series Art. 92). 

If we wish to express F{6) in terms of Zonal Harmonics and to use the 
unit sphere as our region we have 

Ck = f" rF{6) Pi-fcos 6) sin 6 (10 d4> = 27r fF(0) P^icos 6) sin 6 dO. 
A,= J'" JIp,(cos 0)^ sm de del, = 27r[^^J, 

and "* = X ' 

the familiar value. Here 

t^^^ = ^{ 2-/ci^(^)]-^sin ^c^^ -4. 2 2lT-l}- 

(b) Tesseral Harmonics. The Tesseral Harmonics of both tj'pes, i.e., 

cos nif>Pi'i,{iJ.) and sin n<f>P^(fi), a,re a set of functions orthogonal over the 

unit sphere, unless in a pair of the same type the same subscript and tlie same 

index occur. (Byerly's Fourier's Series, Art. 105, and Art. 105, Ex. 2). 

Hence if 

O 



A ' 



where 



C„,^ = f" [fi^' <^) cos n<\> P;(/^) smed0 # , 

A.™ = [J}osH<f>lPMY sin 0d0d^ = ^- [^±i^, 

and b,,^^ = ^, 

where 



6V = fl'/i^^ <!>) si" ^4> ^^(/^) sin d0 d4, ^',,„, = A 



n.m' 



and ao,„ - . 
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where 






then 



m=0 



(v. Byerlj's Fourier's Series, Art. 107). 
Here 



^r 4^ . , ^ 27r (m + n) ! / , ^,>\-|> 
~ 2^ L2m + 1 " -'" "^ 2v 2m + 1 (m - 7i) ! V"'"' ^ "•'«; J ( ' 



15. Miscellaneous Examples. The examples which follow will 
give some notion of the convenience and power of the method of this paper 
in dealing with problems that require the use of the Harmonic Analysis. 

In each of the first three examples a development of unity of importance 
in problems in Mathematical Pliysics is obtained and verified either directly 
or, as in Example 3, incidentally in the solution of an actual Physical Problem. 

The developments in Example 1 and Example 2 are familiar, but are us- 
ually obtained by other methods, that in Example 3 I believe is new. 

In Example 4 an approximate representation of unity in Bessel's Func- 
tions which ai'e not orthogonal is considered. 

In all the examples I refer freely to my Fourier's Series and Spherical 
Harmonics. (Boston, Ginn & Co.) 

16. Example 1. To express 1 as 



S 



S. mirx . n-Kii „ 
^ oin oln y_ * 



a,,, „ sin sin 



Here our functions (^„,_„ (^>2/) are easily seen to be orthocjonal over the 
rectangle bounded by the axes and the lines x = a, y = b. 

* r. Fourier's Series, pp. 127, 128. 
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^ f C" . nirx . niry , , 

^mn = / /sin Sin —5-=^ 0,7/ dx = 0, 

Jo Jo « 



unless m and n are both odd, in which case 

y Aab 

16 



Hence 1 = . 



':;=^';^ l . (2w+l)7ra; . (2w+l)7ry"] 

2v 2^(2m+l)(2«+l) ''" a ^'° 6 • 

L J abl"^-^-:;^ Z^ 2^ {2m+iy(2n + l) 

L m=0 « =0 

but the double series in the last term is easily seen to be the square of the 
series 

111 

12 ^ 32 -1- 52 + 

that is the square of tt^/S, and [JS'^I = and our expression is a correct de- 
velopment of 1 holding good for < ar < a and < 1/ < b. 

17. Example 2. The Bessel's Functions Jo(f^k>')> '^oiw) are orthog- 
onal over the area of the circle whose centre is the origin and whose radius is 
a if and only if fj,f.a and /t^a are roots of the equation Jo(x) = 0, or of 
Ji(x) = 0, or of xJi(x) — \J'o(x) = 0. (Fourier's Series, Art. 125.) 

(a). Let us express 1 as Sa^. Jo{fi.i.r) if fi^a is a root of Jq{x) = 0. 

Here at. = -7^, 

where 

fi" fa ra 27ra 

Cfc = / / Mw) rdrd<f> = 27r / Jo(fif.r)rdr = Ji(Mi.a). 

Jo Jo Jo l^k 

(^Fourier's Series, p. 229 (6).) 
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A, = rr^J,(i.,r)yrdrd<f> = 27rJ^"[jo{ti,r)yrdr= 7ra*[ji(/.,a)Js 



and Ok 






The equation 

J'o(a;) = 1 — -^ + -^^ — 2242^2 + • • • = 

may be regarded as an equation in x^, and the sum of the reciprocals of its 
roots is the negative of the coefficient of its second term. Hence 

and [^*] = ; and (1) is a correct development of 1 over the circle whose 
radius is a and therefore over the range < )• < a. 

(6) If uta is a root of the equation xJ^(^x) — XJi (ic) = 0, we have 
(v. Fourier's Series, p. 229) 



and ar;(. 



2/iij.a Ji(fiii.a) 2\ 






1 - 2\ 'V ''oKf^k') /ON 
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X** 11 



The equation 



. . . = 



may be regarded as an equation in x^. 

The sum of the i-eciprocals of its roots is the negative of the coefficient 
of the second term divided by the first term or (2 + \)/(4\) ; that is, 



2 + X 

l4a^~ 4X ' 



is not so easy. 

We must first transform F(x^) = xJi(x) — \J,)(x) = into an equa 
tion whose roots exceed those of F(x^) = by X*. 

The constant term will be F(— X*) and the coefficient of the term of 
first degree will be F'(— X*). (v. Todhunters' Theory of Equations, §54.). 

Here F'{x^) = ^ [xMx) + XJi(x)] = J [^x) + ?^^] . 

2^x* + / 



Ula^ F{-X^) 2[XtVi(Xi) -XJoCXi)] 



1 Jo(X{) — {Ji(\i) 



2X iJ,(Xt) - Jo{\{) 2X' 



>-A r_ji i__-i ^ 2 + X _ J_ ^ 1 

2~t IfiW 5^* + f4(i^J 4X 2X 4" 
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Hence [^*] = and (2) is a correct development of 1 over the circle 
whose radius is a, and therefore over the i-ango < r < a. 

18. Example 3. A homogeneous cylinder of radius c and altitude b 
is initially at the constant temperature unity throughout. Its entire surface 
is suddenly cooled to the temperature zero and kept there. Required the 
temperature of any point in the cylinder at any time after cooling has begun. 

We have to satisfy the heat conduction equation 

BtU = a^lJD^^u + - I)rU + ^,u\ (1) 

subject to the conditions 

M = when 2 = (2) 

u — Q when 2 = 6 (3) 

t< = when r = c (4) 

M = 1 when ^ = (5) 

It is not difficult to find and it is easy to verify the particular solution 

u = €-(''' + (*') Joifir) sin vz, 

where fi and I'are any constants. This value satisfies equation (1) and con- 
dition (2). 

WITT 

If v = —r- it satisfies (3). If /tt = yitfc, when fi^c is a root of Jo(x) = 0, 
it satisfies (4) . 

Let u =2 2«.,™e -^ '^^''" Mf^^r) sin ^. 

This will be the required solution if 

1 = ]^ ^«*,n. Mf^kr) sin -^ 
throughout the cylinder. 

The functions Jo{t''k>') sin —, — are easily seen to be orthogonal through- 
out the cylinder. 

_ Gk,m 

"le.m — ~3 » 

-^k.m 
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where 

Ci.m = / / / Jo (/**>•) sin ^ r dr d<f> dz 

2cb 
= (1— cos mtr) Ji(fJi'icC) ; (Four. Ser., page 229.) 

A". = ^ ['^i(MfcC)]^ ; {Four. Ser., page 229.) 

4(1— cos m-tr) 

and our expression for 1 is 

4 .r--^ x-\ 1— cos wiTTZ _ , , . mirz 
1 — - > > r^ ^ •^(Mi^) sin — T-. 

r 1 2 1 1 = i :!lf = 

[_p + 32 + 52 + • • "J ^8 



^ (1- cos wtt)" _4rl 24-1-4- 1=i— =!r 

A mV ~ TT I P ''■ 32 "^ 52 ^ ■ ■ J 9r 8 2 



'|(^'--* ('• E.. M«) sn). 



k 

Hence [.E'^] = 0, and our expression for 1 is exact and 

8^"'^ 1 , . (2m + l)7r^ -/t!^'^,A.2, 

rigorously. 

19. Example 4. In the preceding examples the functions in terms of 
which we have had occasion to express unity have in every case been orthog- 
onal over the range or throughout the region considered. As an example of 
a different sort suppose we wish to express 1 in terms of Bessel's Functions 
J^^jjkjpc) which are not orthogonal over the circle over which the expression is 
to be used.* 

* For a case of this sort see a paper by B. O. Peirce, Proc. Amer. Acad., vol. 43, No. 6, 
Sept. 1907, p. 175, pp. 180-184. 
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There is no serious diffioulty in getting any specified number of terms 
by the formulas of §6 of this paper. 
We have 

Oi = / Jo{nir)r dr d^ = — Ji(f^ia) ; 

Jo J f'l 

A =fJJ LMt'tr)]' rdrd<l> = ira^ \ [Jo(l*i<i)f + lJi(^iia)f\ ; 

/2t fa 

f''k~ri 

(v. Fourier's Series, §§125-126, pp. 228-229.) 
For instancCj if a single term ctiJoif^ir) will suffice, 



Mia U JiCMi")) + \Ji{M-ia)j J 



r^*1 = 1 - 



4[Ji(/^ia)]^ 



Mi«' [(^o(Mia)) + (jiCiCtia))'] 

If we need three terms we have 

AiUi + B^^^a^ + i?i,3a3 = d, 

J^lfi^l "t" -"2^2 "H -"2,3% ^ ^2> 
■^1,3*1 + B<i^i + -^303 = Cg, 

and 

L^ J = :^ I 'r«^ - -4i«i - -djol - A^al — 2 r2?i_2aia2 + Bi/t^a^ + -^Lsaifls 1 • 

20. Conclusion. It may be objected that the test of excellence used 
in this paper is a purely artificial one. It is, however, often a very useful 
one and for many purposes it is precisely what is wanted ; and it is in fact 
the test that a specified number of terms in any of the received Harmonic de- 
velopments always fulfil. 
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The method given here like other rough and ready methods, for instance 
that of Lagrange (Fourier^s Series, §§19-25, et passim), has nothing to do, 
except indirectly, with the question of development in series. Indeed if a 
development in series of the proposed form is possible the coefficients of a 
specified pet of terms given by the method of Least Squares are not the coeffi- 
cients of the corresponding terms of the series unless the functions happen 
to be orthogonal throughout the region considered. In other cases they 
give a closer approximation than would the actual terms of the series. 
(V. Ex. 1 (c) §6.) 

The check on the accuracy of the representation in any concrete case 
given by the value of [-E"^] is useful even when terms of a true series devel- 
opment are employed and it is easily applied in any case since all the con- 
stants in [^^] must have been calculated already in getting the coefficients. 

If the functions used in the representation of f(x) are orthogonal 
throughout the region in question the labor of getting the coefficients is much 
lessened, since every one is obtained separately, and every additional one ob- 
tained adds to the accuracy of the result. 

If the functions employed are not orthogonal throughout the region the 
coefficients obtained are for many purposes the best possible, and whether or 
not the representation is accurate enough for the matter in hand is at once 
shown by the value of [^*] . 

In the subjects which gave rise to the Harmonic Analysis, i.e., in the 
Motions of Elastic Strings and in the Conduction of Heat, most of the devel- 
opments required are in terms of functions orthogonal over the range or 
throughout the region concerned, and the old theory worked well enough. 
In problems in Electricity and Magnetism, and in Hydromechanics similar 
developments where the functions are not orthogonal are often desired and to 
them the old theory gives no clue. 

Although our method has nothing to do with development in series ex- 
cept indirectly, corollaries (b) and (c) §11 show that it may often be of great 
service in that connection, and they suggest some beautiful little problems 
for the workers in the science of Integral Equations. 

Harvard University, 
Cambridge, Mass., 
July, 1910. 



